Single crystal growth and physical properties of the layered arsenide BaRh2As2 
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Single crystals of BaRli2As2 have been synthesized from a Pb flux. We present the room temper- 
ature crystal structure, single crystal x-ray diffraction measurements as a function of temperature 
T, anisotropic magnetic susceptibility x versus T, electrical resistivity in the a6-plane p versus T, 
Hall coefficient versus T and magnetic field H, and heat capacity C versus T measurements on the 
crystals. The single crystal structure determination confirms that BaRh2As2 forms in the tetragonal 
ThCr2Si2 type structure (space group I4/mmm) with lattice parameters a = h = 4.0564(6) A and 
c = 12.797(4) A. Band structure calculations show that BaRh2As2 should be metallic with a small 
density of states at the Fermi energy N{Ef) — 3.49 states/eV f.u. (where f.u.= formula unit) for 
both spin directions. p{T) data in the a&-plane confirm that the material is indeed metallic with a 
residual resistivity p(2 K) = 29 /xSl cm, and with a residual resistivity ratio p(310 K)/p(2 K) = 5.3. 
The observed x(^) is small (~ 10~® cm^/mol) and weakly anisotropic with Xab/Xc ~ 2. The C{T) 
data indicate a small density of states at the Fermi energy with the low temperature Sommerfeld 
coeflicient 7 — 4.7(9) mj/mol K^. There are no indications of superconductivity, spin density wave, 
or structural transitions between 2 K and 300 K. We compare the calculated density of states versus 
energy of BaRh2As2 with that of BaFe2As2. 



I. INTRODUCTION 

The recent discovery of high temperature supercon- 
ductivity in LaFeAs0i_2;Fa; with Tc = 26 K (Ref. Ij) has 
generated a lot of activity in the search for related super- 
conductors. Even higher TcS of 35 K to 55 K have been 
achieved when La is replaced by other rare-earth elements 
R = Ce, Pr, Nd, Sm, Gd, Tb, and Dy.l^^^^^^^ These 
materials form in the tetragonal ZrCuSiAs type structure 
{P4/nmm)^ The structure is made up of FeAs layers al- 
ternated by LaO layers stacked along the c-axis. The 
undoped materials i?FeAsO show a spin density wave 
(SDW) transition and a coupled structural transition at 
high temperatures (T ^150 K)i^ii£iii On doping with F, 
the SDW and the structural transitions are suppressed 
and superconductivity is observed i^i^i^i^iSiiiSiiiii^ 

Recently a different family of compounds y4Fe2As2 
{A = Ba, Sr, Ca, and Eu) was discovered which crys- 
tallize in the tetragonal ThCr2Si2 type structure and 
has similar FeAs layers seperated by A layers stacked 
along the c-axis. Like the i?FeAsO materials, these 
materials also show SDW and structural transitions at 
high temperatures. ^^ii^^iii^iia^^^ When the A 
atoms are partially replaced by K, Na, or Cs, the SDW 
and structural transitions are suppressed and supercon- 
ductivity is observed The common feature in 
all the materials mentioned above is the FeAs layer in 
their crystal structure. It is of interest to look for other 
materials with related structures and investigate their 
physical properties to see if these can be potential parent 
compounds for new high temperature superconductors. 

The series of compounds ARh2As2 {A = Ba, Sr, and 
Eu) were previously synthesized in single crystal form 
and their crystal structure was reportedi^ The Ba and 
Eu compounds form in the ThCr2Si2 type structure while 



the Sr compound forms in a different structure. The 
crystal structure of BaRh2As2 is shown in Fig. [TJ The 
structure is built up of RhAs layers alternated by Ba lay- 
ers stacked along the crystallographic c-axis. Within the 
RhAs layers the As atoms lie outside the plane formed by 
the Rh atoms. To the best of our knowledge the physical 
properties of these v4Rh2As2 compounds have not been 
investigated before. 

Herein we report on the crystal growth, single crystal 
structure, resistivity p in the a5-plane versus temperature 
T, Hall effect versus T and magnetic field H, magnetic 
susceptibility x versus T, and heat capacity C versus T 
measurements of BaRh2As2. Our experimental results 
are compared with predictions of band structure calcula- 
tions. 



II. EXPERIMENTAL DETAILS 

Single crystals of BaRh2As2 were grown out of 
Pb flux.^^ The elements were taken in the ratio 
Ba:Rh:As:Pb = 1.1 : 2 : 2.1 : 50, placed in an alu- 
mina crucible and then sealed in a quartz tube under 
vacuum (10~^ mbar). The whole assembly was placed 
in a box furnace and heated to 1000 °C at a rate of 
50 °C/hr, left there for 10 hrs and then cooled to 500 °C 
at a rate of 5 °C/hr. At this temperature the molten Pb 
flux was decanted. Shiny plate- like crystals of typical size 
1.5 mmxl.5 mmxO.l mm were obtained. The composi- 
tion of the crystal was checked using energy dispersive x- 
ray (EDX) semiquantitative analysis using a JEOL scan- 
ning electron microscope (SEM). The SEM scans were 
taken on a cleaved surface of a crystal. The EDX gave 
the average elemental ratio Ba:Rh:As = 20.4 : 40.8 : 38.8 
which is consistent with an approximate 1:2:2 stoichiom- 
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FIG. 1: The crystal structure of BaRh2As2. Layers of Ba 
and layers of RhAs alternate along the c-axis. Within the 
RhAs layers the As atoms lie outside the plane formed by the 
square-lattice Rh atoms. 



etry for the compound. An SEM image of a typical as- 
grown crystal is shown in Fig. ^a). Laue back scat- 
tering measurements on the crystals showed that the 
largest surface of the plates was perpendicular to the c- 
axis. For crystal structure determination, single crystal 
x-ray diffraction measurements were done at temperature 
T = 293 K using a Bruker CCD-1000 diffractometer with 
Mo Ka (A = 0.71073 A) radiation. To look for struc- 
tural phase transitions, single crystal x-ray diffraction 
measurements were performed on a standard four-circle 
diffractometer using Cu Ka radiation from a rotating an- 
ode x-ray source, selected by a Ge(l 1 1) monochroma- 
tor. For these measurements, a plate-like single crystal 
with dimensions of 2x1.5x0.2 mm^ was selected. The 
sample was mounted on a flat copper sample holder in 
a closed cycle displex cryogenic refrigerator with the (0 
1)-(1 1 0) reciprocal lattice plane coincident with the 
scattering plane. The diffraction patterns were recorded 
for temperatures between 10 K and 300 K and with a 
setup optimized for high resolution in transverse scans 
in the scattering plane. The measured mosaicity of this 
crystal was 0.025° full width half maximum for both the 
(0 10) and (1 1 10) reflections at room temperature 
as shown in Fig. ^h) , indicating the good quality of the 
single crystal. 

The x(T) was measured on a collection of aligned crys- 
tals with a total mass of 5.11 mg using a commercial 
Superconducting Quantum Interference Device (SQUID) 
magnetometer (MPMS5, Quantum Design), the standard 
four-probe p{T) was measured with a current of am- 
plitude 1 = 2 mA at a frequency of 16 Hz, using the 
ac transport option of a commercial Physical Property 
Measurement System (PPMS5, Quantum Design). The 
contacts were made with silver epoxy on a cleaved sur- 
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FIG. 2: (a) A scanning electron microscope (SEM) image of 
a typical as grown crystal of BaRh2As2. The white mottled 
area on the surface is a thin layer of Pb from the flux, (b) 
Rocking scan through the (0 10) reflection, and (c) rocking 
scan through the (1 1 10) reflection of a BaRh2As2 single 
crystal. 



face of a crystal. The current 1 = 2 mA was applied in 
the a6-plane. The C{T) was measured on a collection of 
three crystals of total mass 3.2 mg using the commercial 
PPMS. 

For the density of states (DOS) calculations, we have 
used the full potential linearized augmented plane wave 
(FP-LAPW) method with a local density approxima- 
tion (LDA) functionali^ The difference in energy of 
0.01 mRy/cell between successive iterations was used as a 
convergence criterion. The value of -RMT-K^max [the small- 
est muffin tin (MT) radius multiplied by the maximum 
k value in the expansion of plane waves in the basis set] 
which determines the accuracy of the basis set used, is 
set to 9.0. The total Brillouin zone was sampled with 405 
fc-points in the irreducible Brillouin zone. The employed 
MT radii are 2.4, 2.2 and 2.2 atomic units for Ba, Rh, and 
As, respectively. The structural data (lattice constants 
and atomic positions) were taken from the previously re- 
ported single crystal structure of BaRh2As2'2i 
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III. RESULTS 

A. Single Crystal X-ray diffraction and Structure 
of BaRh2As2 

A small (0.35 mmxO.35 mmxO.07 mm) plate-like sin- 
gle crystal was used for crystal structure determination. 
The initial cell parameters were taken as those previ- 
ously reported for BaRh2As2 (ThCr2Si2 structure, Z = 2 
formula units/unit cell, space group IA/m'mm)M' The 
final cell parameters and atomic positions were calcu- 
lated from a set of 861 strong reflections with good 
profiles in the range 29 = 6°-61°. The unit cell pa- 
rameters were found to be a = 6 = 4.0564(6) A and 
c = 12.797(4) A. These values are in good agreement 
with previously reported values for single crystalline 
BaRh2As2 [a = 6 = 4.053(1) A, and c = 12.770(3) A]^ 
There is only one atomic coordinate not constrained by 
symmetry requirements, the z position for As. Our 
value of z = 0.3566(3) is also in good agreement with 
z = 0.3569(1) reported previously by the full single crys- 
tal refinement. The single crystals of BaRh2As2 have a 
high tendency to split/cleave into very thin plates per- 
pendicular to the c-axis resulting in significant mosaic- 
ity/twinning. This resulted in some reflections to be 
broad/extended. This and the unfavorable aspect ratio 
of the crystal for the appropriate absorption correction 
did not allow a full refinement of the structure with a 
reasonable R factor. 

Figures ^a,) and (b) show the results of our tempera- 
ture dependent single crystal x-ray diffraction measure- 
ments which were done to check for the possibility of 
a structural phase transition in BaRh2As2. Figure [3]^a) 
shows the longitudinal (0 scans through the (0 
10) reflection at various temperatures between 10 K and 
300 K relative to the ahgnment at T = 200 K. We did not 
observe any change in the shape of the (0 10) reflection 
as function of temperature. The change in the position of 
the peak results mainly from lattice parameter changes 
between different temperatures. The small peak at the 
shoulder of the (0 10) reflection at low temperatures 
[Fig. El^a)] arises from a second grain with a slightly dif- 
ferent orientation. In contrast to the isostructural com- 
pounds AFe2As2 [A = Ba, Sr, Ca,) ^i^S,'-^'' transverse (^ 
^ 0) scans through the (1 1 10) reflections [Fig. [^b)] 
show no splitting and/or abrupt changes in reflection po- 
sition over the temperature range between 10 and 300 
K, indicating the absence of the structural transition re- 
ported for the isostructural AFe2As2 {A = Ba, Sr, Ca) 
compoundsji^ii^iii This is consistent with the macro- 
scopic measurements discussed later. 



B. Density Of States 

The calculated total density of states (DOS) N{Ef) 
for BaRh2As2 for both spin directions versus energy E 
measured relative to the Fermi energy E-p, is shown in 
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FIG. 3: (a) Longitudinal (0 Q scans through the position 
of the (0 10) reflection for selected temperatures. The off- 
set between each data set is 1000 counts/s. (b) Transverse 
{£, C 0) scans through the position of the (11 10) reflection 
for selected temperatures. The offset between each data set 
is 500 counts/s. In (a) and (b), the horizontal axes represent 
the scattering vectors relative to the alignment at T = 200 K 
to illustrate the temperature dependence of the reflection po- 
sitions. 



Fig. IHa). We calculate N{Ef) = 3.49 states/eV f.u. 
(f.u. = formula unit) for both spin directions for 
BaRh2As2. Figure [DJb) shows the partial DOS for the 
individual atoms. The maximum contribution to the to- 
tal N{Ef) is from the Rh 4c? states with small but signif- 
icant contributions from Ba and As. Since rather small 
MT radii were employed, wave functions for the 4c? Rh 
states extend out into the interstitial region. This results 
in a significant DOS in the interstitial region which is not 
accounted for in the partial DOS estimates. This inter- 
stitial DOS is however accounted for in calculation of the 
total DOS. Therefore, the partial DOS shown in Fig.lHb) 
do not add up to the total DOS shown in Fig. Ufa). 
It is of interest to compare the DOS versus energy for 
BaRh2As2 with that of BaFe2As2. Figure [5] shows the 
total DOS for both spin directions for BaRh2As2 (solid 
curve) and BaFe2As2 (dashed curve). The band struc- 
tures for the two materials are seen to be qualitatively 
quite similar. It appears that the bands for BaRh2As2 
are shifted down (and stretched) in energy relative to the 
bands for BaFe2As2. For BaFe2As2, if the Fermi energy 
Ep moves down by about 0.5 eV on hole doping, there 
is a large increase in the DOS (almost 4-5 times) com- 
pared to the parent compound. To move the E-p for the 
BaRh2As2 material to the corresponding band located at 
about —2 eV below Ep would evidently require a much 
larger amount of doping. 
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FIG. 4: (a) The total density of states for both spin directions 
for BaRh2As2. (b) Partial density of states of Ba, Rh, and 
As atoms in BaRh2As2. 



C. Magnetic Susceptibility 

The magnetic susceptibilities for BaRh2As2 measured 
with an applied magnetic field H ~ 2 T parallel to 
the c-axis Xc{T) and with H parallel to the a6-plane 
XahiT) are shown in Fig. [6] (after correcting for con- 
tributions from the sample holder). The measurements 
were made on a collection of four c-axis aligned crys- 
tals of total mass 5.11 mg. Both Xc{T) and Xab(T) 
are almost temperature independent and have very 
small average values Xc{T) = 1.5xl0~^ cm'^/mol and 
Xah{T) = 3xl0~^ cm'^/mol. There is a small anisotropy 
in the whole temperature range with Xab/Xc ~ 2. The 
data are somewhat noisy, most likely due to both the 
small absolute values of the susceptibility and the small 
crystal mass used for the measurements. The powder av- 
eraged susceptibility x{T) can be calculated for a tetrag- 
onal system as 



X 



2Xab + Xc 



(1) 



From the average values of Xab and Xc we estimate an 
average value x = 2.5x 10~^ cm'^/mol from 10 K to 300 K. 
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FIG. 5: Comparison of the total density of states for both spin 
directions for BaRh2As2 (solid curve) with those for BaFe2 As2 
(dashed curve). 



The susceptibility x{T) can be written as 



X = Xcoro + XL + XVV + XP 



(2) 



where Xcorc is the diamagnetic orbital contribution from 
the localized core electrons (ionic or atomic), xl is the 
Landau orbital diamagnetism of the conduction elec- 
trons, XVV is the Van Vleck paramagnetic orbital con- 
tribution and xp is the Pauli paramagnetic spin suscep- 
tibility of the conduction electrons. We evaluated xp and 
Xvv from electronic structure calculations using^ 



XP = AN{Ey) 



(3) 



where /ie is the Bohr magneton and N{E-p) is the bare 
density of states at the Fermi level for both spin direc- 
tions, and 
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(4) 

where / is the Fermi-Dirac distribution function and Lz 
is the z component of the angular momentum operator 
We have calculated the anisotropic xvv for magnetic field 
H in the (001) and (110) directions. 

The diamagnetic contribution xd — Xcoro + Xl to 
the total susceptibility was also estimated. Since dia- 
magnetic susceptibility is problematic to calculate in ex- 
tended systems because of inter-cell currents, we take the 
value from calculations of the neutral atoms and ignore 
the solid state effects. This gives a lower limit estimate of 
Xd- It also misses any anisotropy arising from inter-cell 
currents. 

The results of these calculations are summarized in Ta- 
ble |T1 It can be seen that although the individual contri- 
butions Xp , Xd , and xvv are large there is a large cancel- 
lation between the positive xp and xvv , and the diamag- 
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TABLE I: Estimated contributions to the total susceptibility 
X in units of 10~^ cm^/mole 



contribution 


Value 


XP 


112 


Xd 


-222 


Xvv 


H II (001) H II (110) 
153 138 




300 



FIG. 6: The magnetic susceptibility x versus temperature T 
for BaRh2As2 for the magnetic field H applied parallel and 
perpendicular to the tetragonal c-axis. 
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FIG. 7: The resistivity p versus temperature T in the afe-plane 
for a crystal of BaRh2As2 between 2 K and 300 K. The inset 
shows the p{T) data between 2 K and 22 K on an expanded 
scale to highlight the residual resistivity po and the slight 
upturn below about 5 K. As noted, the residual resistivity 
ratio is RRR = /9(310 K)/p{2 K) = 5.3. 



is intrinsic to the material. The residual resistivity ra- 
tio RRR = /9(310 K)/p{2 K) = 5.3 for our BaRhsAss 
crystal is typical of what is observed for single crystals 
of the isostructural y4Fe2As2 {A = Ba, Sr, Ca, and Eu) 
materials^^>i^^>28^ 



netic xd which results in a small net x- Although the to- 
tal X from the calculations is of the same order as the ex- 
perimentally measured value, the calculations could not 
quantitatively reproduce the experimental numbers and 
the anisotropy (Fig. [6]). Since these calculations did not 
take into account the anisotropy of xd associated with 
valence electrons, the above comparison strongly suggests 
that xd also has a large anisotropy, with a larger value 
for fields along the c-axis. 



D. Resistivity 

The electrical resistivity p versus temperature T in 
the a6-plane is shown in Fig. [71 The p{T) data show 
typical metallic behavior with a monotonic decrease 
with decreasing T between 310 K and 10 K. Below 
T = 10 K p{T) saturates to a residual resistivity value 
p(2 K) = 29(3)/Ltr2 cm. There is however a slight up- 
turn in p{T) at the lowest temperature as seen in the 
inset in Fig. [71 A corresponding upturn is also observed 
at low temperatures in the heat capacity C{T) data dis- 
cussed below. This upturn in both the p{T) and C{T) 
data at low temperatures may be due to the presence of 
a small amount of paramagnetic impurity in the mate- 
rial. However, it cannot be ruled out that this behavior 



E. Hall effect 

Figure [SJa) shows the dependence of the Hall coeffi- 
cient i?H on magnetic field H at various temperatures. 
The i?H is negative at all fields and temperatures which 
indicates that electron conduction dominates the elec- 
tronic transport. The i?H is approximately constant with 
magnetic field at all temperatures. From i?H one obtains 
a lower limit estimate of the carrier density using the 
single-band expression^ 



where n is the carrier density and q is the charge of the 
current carriers. Using the value i?H = 3.5 x 10~^ cm'^C"^ 
at 300 K and 8 T and q = — e where e is the magnitude 
of the charge of the electron, we get n = 1.8x10^^ cm"'^. 
One can get an estimate of the density of states at the 
Fermi energy N{E-p) using the single-band relation^ 

N{Ey) = — , (5) 

where m is the free electron mass, Vf.u. = 210.5 A"^ is the 
volume per formula unit, and h is Planck's constant di- 
vided by 27r. With the value of n obtained above one gets 
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mass m = 3.2 mg of the crystals used for the measure- 
ment. However, the heat capacity at room temperature 
C(300 K) « 120 J/mol K is close to the classical Dulong 
Petit value C = 125 J/mol K expected for BaRh2As2. 
Below T — 100 K the quality of the data is better. There 
is no signature of any phase transition in the temperature 
range of our measurements. 

FigurelUb) shows the C{T)/T versus data between 
2 K and 20 K. The C{T)/T data show curvature in the 
whole temperature range shown. This suggests that there 
may be anharmonic contributions to the lattice heat ca- 
pacity even at low temperatures. Only below about 10 K 
do the C{T)/T data show a quasi-linear dependence on 
T^. At the lowest temperature the C /T shows a tendency 
of saturating to a value C(2 K)/T = 15.5 mJ/mol 
as shown in the inset in Fig. [IJb) where the C/T ver- 
sus T data are shown below T = 6 K. From fits to 
the C/T versus T data up to 10 K by the expression 
C/T = j+PT'^, where 7 is the Sommerfeld coefficient of 
the electronic heat capacity and /3 is the coefficient of the 
lattice heat capacity, we obtain 7 — 4.7(9) mJ/mol 
and P = 1.93(4) mJ/mol K^. 

The density of states at the Fermi energy for both spin 
directions N{Ep) can be estimated from the above value 
of 7 using the relation^ 



FIG. 8: (a) The Hall coefficient Ru versus magnetic field H 
measured at various temperatures, (b) The Hall coefficient 
Ru versus temperature T in a magnetic field of 8 T. 



N{Ep) = 2.3 states/(eV f.u.) for both spin directions 
which is smaller than the 3.49 states/(eV f.u.) obtained 
above from band structure calculations. This suggests 
that both electrons and holes contribute to the electrical 
conduction and the net carrier density estimated from 
the Hall measurements is smaller than the actual value 
for the material, and/or that the conduction electron ef- 
fective mass is greater than the free electron mass. 

Figure [5]Jb) shows the plot of i?H versus temperature 
T at a magnetic field of 8 T. The i?H is almost con- 
stant between T = 50 K and 300 K. Below T = 50 K 
however, it decreases strongly by a factor of 4-5. In a 
single band model the Ru is expected to be temperature 
independent if the scattering rate r is isotropic. How- 
ever, in a two-band model, the i?H could be temperature 
dependent if, for example, the scattering rates of the two 
bands have a different temperature dependence or the 
fractions of charge carriers in the two bands change with 



(6) 



where Acp is the electron-phonon coupling constant. 
Setting Acp = as a first approximation, one gets 
N{Ef) = 2.0(4) states/(eV f.u.). This is a smaller 
than the value predicted from band structure calcula- 
tions N{Ef) = 3.49 states/(eV f.u.). 

From the value of /3 = 1.93(4) mJ/mol K"* estimated 
above one can obtain the Debye temperature Od using 
the expression^ 



5P 



1/3 



(7) 



temperature 



33 



where R is the molar gas constant and n is the number 
of atoms per formula unit {n — 5 for BaRh2As2). We 
obtain Od = 171(2) K for BaRh2As2. The values of 7 
and 8d estimated for BaRh2As2 are similar to what has 
been reported for the isostructural compounds AFe2As2 
{A = Ca, Sr, Ba, and Eu)^i^^ii^ 



IV. CONCLUSION 



F. Heat Capacity 

The heat capacity C versus temperature T of 
BaRh2As2 between 1.8 K and 300 K is shown for a col- 
lection of three crystals in Fig. [DJa). The data at high 
temperatures are somewhat noisy because of the small 



We have synthesized single crystalline samples of the 
layered rhodium arsenide BaRh2As2 and characterized 
them using single crystal x-ray diffraction, anisotropic 
magnetic susceptibility versus temperature, resistivity 
versus T, and heat capacity versus T measurements be- 
tween 2 K and 300 K. The single crystal structure de- 
termination confirms that BaRh2As2 crystallizes in the 



7 




50 100 150 200 250 300 
T(K) 



tertagonal ThCr2Si2 type structure with lattice parame- 
ters a = b = 4.0564(6) A and c = 12.797(4) A. Single 
crystal x-ray diffraction meaasurements down to 10 K did 
not show any evidence for a structural transition. The 
density of states calculations give a total density of states 
for both spin directions N{Ey) = 3.49 statcs/(cV f.u.) 
with maximum contribution from the Rh Ad states. The 
x(T) is small and temperature dependent. The small 
x(T) suggests a small density of states. This is supported 
by the C(T) data which give a small Sommerfeld coeffi- 
cient 7 = 4.7(9) mJ/mol K^. 
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FIG. 9: (a) Heat capacity C versus temperature T for 
BaRh2As2 crystals between 2 K and 300 K. (b) C/T versus 
data between 2 K and 20 K. The inset shows the C{T)/T 
versus data below T = 6 K. The solid line through the 
data in the inset is a fit by the expression C/T = j+^T^. 
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